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We calculate the amplitudes involved in the heavy baryon nonleptonic decay Kh KJ/'tp using 
perturbative QCD factorization theorem, which are expressed as convolutions of hard b quark decay 
amplitudes with the Kb baryon, A baryon and J/ip meson distribution amplitudes. It is found that 
nonfactorizable contributions dominate over factorizable ones. Because of soft cancellation in pairs 
of nonfactorizable diagrams, the Kb KJ/ij) decay is characterized by a large scale, such that 
perturbative QCD is applicable. Employing the distribution amplitudes determined in our previous 
works and from QCD sum rules, we derive the branching ratio B{Kb —> KJ/ip) = (1.7 ~ 5.3) x 10“"* 
in agreement with data. We predict an asymmetry parameter a = —0.17 ~ —0.14 associated with 
the anisotropic angular distribution of the A baryons produced in polarized Kb baryon decays. 


I. INTRODUCTION 

Recently, we have developed perturbative QCD (PQCD) factorization theorem for the semileptonic heavy baryon 
decays Kb plv and A;, —> Kclv [^,H. This theorem states that nonperturbative dynamics involved in a high- 
energy QCD process can be factorized into hadron distribution amplitudes, and the remaining piece is calculable in 
perturbation theory. A distribution amplitude, though not calculable, is universal. Once a distribution amplitude 
is determined from experimental data of some processes, it can be employed to make predictions for other processes 
involving the same hadron. According to this theorem, the form factors involved in the decay Kb —> KJ,i> have been 
expressed as the convolutions of hard b quark decay amplitudes with the universal Kb and A^ baryon distribution 
amplitudes. It has been found that perturbative contributions to the Kb —> decays become more important at the 
maximal recoil of the Ac baryon with the velocity transfer about 1.4. This observation indicates that PQCD is an 
appropriate tool for analyses of two-body nonleptonic Af, baryon decays. 

An essential feature of PQCD is that it goes beyond the conventional approach to exclusive nonleptonic heavy 
hadron decays based on the factorization approximation (FA) . Both factorizable and nonfactorizable contributions 
from various topologies (emission and annihilation) can be evaluated systematically in PQCD. Though nonfactorizable 
contributions are usually negligible, they become dominant in the modes B —> J/ipK^ *) @ , whose factorizable 
contributions arise from internal IF-emission with the small Wilson coefficient 02 defined below. This is the reason 
it is difficult to accommodate the data of the ratios R = B{B J/tpK*)/B{B J/ijjK) and Rl = B{B 

J/il^K'l)/B{B J/tpK*) simultaneously in FA. However, these data can be explained in the PQCD formalism Q|. 

In this paper we shall extend PQCD factorization theorem to the more complicated baryon decay Af, ^ AJ/'0, and 
show that nonfactorizable contributions also play an essential role. A simple investigation indicates that each diagram 
for the Af, ^ AJ/ijj decay is characterized by the scale A, where A = Ma^ — mf, is the mass difference between the 
Af, baryon and the b quark. Factorizable contributions, having this low typical scale, are suppressed by the Wilson 
coefficient 02 . The dominant nonfactorizable contributions, due to soft cancellation between a pair of diagrams, are 

characterized by the higher scale ^AMa^[1 — with Mj/^ being the J/ip meson mass. Therefore, the 

Af, —> KJjil) decay is a special heavy-to-heavy mode, to which PQCD is applicable. 

The Af, baryon distribution amplitude has been chosen to satisfy the experimental upper bound of the At, —> Kclv 
branching ratio and the requirement of heavy quark symmetry 0. The J/^/> meson distribution amplitudes have 
been determined from the experimental data of the B J decays 0. The A baryon distribution amplitudes 
have been derived from QCD sum rules |0,0. We shall employ these distribution amplitudes, because of their 
universality, to predict the branching ratio B\Kb AJ/ijj). We also consider another interesting quantity, the 
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asymmetry parameter a defined in Sec. V, which is related to the anisotropic angular distribution of the A baryons 
produced in polarized A;, baryon decays. 

The Af, —> AJ/t/j decay has been discussed in the quark model based on FA Q. Our predictions 

B{Ab AJ/^) = (1.7 ~ 5.3) X 10"^ , 

a =-0.17 ~-0.14 , (1) 

are consistent with those from FA. However, the theoretical bases of the two approaches are quite different: nonfac- 
torizable contributions are neglected, and the Wilson coefficient 02 is treated as a free parameter in FA. We emphasize 
that there is no free parameter in the PQCD calculation. The comparison of our predictions with data will provide a 
justification of the PQCD formalism for heavy baryon decays. 

In Sec. II we briefly review factorization theorem for exclusive nonleptonic heavy hadron decays. The hadron 
distribution amplitudes and the associated Sudakov factors are defined in Sec. III. The factorization formulas for 
the Ab AJ/tp decay amplitudes are presented in Sec. IV. Numerical results are discussed in Sec. V. Section VI is 
the conclusion. To deliver our ideas clearly, we shall leave all technical calculations and complicated expressions to 
the Appendices. Appendix A contains the useful integrals employed in this work and the measures of phase-space 
integrations for different diagrams. The hard amplitudes are summarized in Appendices B and C. 


II. FACTORIZATION THEOREM 


We briefly review PQCD factorization theorem for exclusive nonleptonic heavy hadron decays Q, and discuss its 
application to the A;, ^ AJ/ip mode. The effective Hamiltonian responsible for this mode is expressed as 


+ C2(m)02(m)] , 


( 2 ) 


where Gp is the Fermi coupling constant, V’s the Cabibbo-Kobayashi-Maskawa matrix elements, Cl^ 2 {^J■) the Wilson 
coefficients, and fi an arbitrary renormalization scale. The four-fermion operators Oi ^2 are written as 


Oi = (cb)(sc) , O 2 = (sb)(cc) , 


( 3 ) 


with ( 9192 ) = 9i7m(1 - 75)92 being the V - A current. 

We shall have a careful look at the derivation of the above effective Hamiltonian starting with eight-quark amplitudes 
(three quarks from the A;, baryon, three quarks from the A baryon and two quarks from the J/ip meson). The 
lowest-order diagrams contain one W boson for the weak decay of the b quark and two hard gluons attaching the 
spectator quarks. Nonleptonic A& baryon decays involve three scales: the W boson mass Mw, at which the matching 
conditions of the effective Hamiltonian are defined, the hard scale t related to , which reflects specific dynamics of 
different modes, and the factorization scale of 0(A) Q. The factorization scale is introduced to separate perturbative 
contributions from nonperturbative contributions absorbed into hadron distribution amplitudes (p(x,b), x being a 
momentum fraction carried by one of the valence quarks. In the PQCD formalism the factorization scale is chosen as 
I/&, b being the transverse extent of a hadron. 

Radiative corrections generate various types of large logarithms, such as as ln(Mw/t) and Ug ln(t6), which should 
be summed by renormalization-group (RG) methods to give evolution factors. The Wilson coefficients C{t) from the 
summation of the first type of logarithms correspond to the evolution from Mw to t. The second type of logarithms 
is summed to give the evolution g{t,b) from t to I/&, which is governed by an anomalous dimension different from 
that of C{t). The difference arises from the loop corrections to the eight-quark amplitudes for g{t,b), and those to 
the four-quark amplitudes associated with the W boson exchange for C(t) 0 - There also exist double logarithms 
Us ln^(P6), P being the dominant light-cone component of hadron momentum, which appear in radiative corrections 
to hadron wave functions. These logarithms, from the overlap of collinear and soft divergences, are treated by the 
resummation technique |I5 The result is a Sudakov exponential exp[—s(P, &)], which decreases fast with b and 


vanishes at 6 = I/Apcp Aqc d being the QCD scale. Since the Sudakov factor supresses long-distance contributions 
from the large b region [p^ , |l8| , the hard scale t, always larger than the factorization scale 1/b, does not go down to A. 

After summing logarithmic corrections, the hard amplitude H(t) can be calculated perturbatively by means of 
Feynman diagrams with on-shell external quarks. The on-shellness in the present leading-power analysis means that 
the virtualities of the external quarks are at most of 0{A/mb). The A;, ^ AJ/ip hard amplitude contains both 
factorizable and nonfactorizable internal kF-emission contributions from the lowest-order diagrams in Fig. I. Note 
that we do not exhibit the diagrams, which are equivalent under exchange of the u and d quarks. It will be observed 
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that nonfactorizable diagrams, especially those with both the charm quarks in the J/il> meson attached by the hard 
gluons, dominate. The vertex corrections to the four-fermion operators, which are calculated in the improved QCD 
factorization [^9[j2l], are of higher orders in the PQCD approach and not considered in the present leading-order 
formalism. 

A naive power counting for the A;, —> KJ/ip decay indicates that the important kinematic region corresponds to the 
parton longitudinal and transverse momenta of 0(A). This is the reason each diagram for this mode is characterized 
by the scale A as stated in the Introduction. The factorizable contributions, characterized by this low scale, are 
suppressed by the small Wilson coefficient a 2 = C 2 + Ci/Nc, Nc being the number of colors. The Ab AJ/ip 
decay is then dominated by the nonfactorizable contributions as in the B J/'ipK^*'> decays. Since the J/ip meson 
is a color-singlet object with the two valence quarks moving parallelly, there exists soft cancellation in a pair of 
nonfactorizable diagrams such as Figs. 1(f) and 1(h). Hence, the parton momenta in the A baryon become of 

0{Mj^^ [1 — and the A{, ^ AJ/ip decay is characterized by the larger scale 

It then makes sense to apply PQCD to this heavy-to-heavy decay mode. 

At last, the Ab —*■ AJjpj decay amplitudes are expressed as the convolution of the above factors, 

C{t) ^ H{t) ^ g{t, b) ® exp[—s(P, b)] 0 ^{x, b) . (4) 

All the convolution factors, except 4>{x,b), are calculable. The hadron distribution amplitude (j){x,b), though not 
calculable, are universal, since they absorb long-distance dynamics of a decay process, which is insensitive to short- 
distance dynamics involved in the b quark decays. Based on universality, we employ the hadron distribution amplitudes 
extracted from other experimental data or from QCD sum rules to make predictions for the A;, —> AJ/ip decay. Note 
that the hard scale t is a convolution variable, since the derivation of Eq. (^ starts with the eight-quark, instead of 
four-quark, amplitudes. Its precise expression should be determined by diminishing next-to-leading-order corrections 
to the hard amplitude H{t). In this work we shall consider the different choices of t as the main source of the 
theoretical uncertainty. 


III. DISTRIBUTION AMPLITUDES 


We define kinematics of the initial and final hadrons as follows. The A^ baryon is assumed to be at rest, and the 
A baryon, regarded as being massless, recoils in the minus direction. The momenta p, p' and q = p — p' of the A;, 
baryon, the A baryon and the J/tp meson, respectively, and the momenta of their valence quarks are parametrized as 

p= (p+,p-,0) = ^(1,1,0) , 

fci = (a;ip+,p“,kiT) , fe = (a;2P^,0,k2T) , fca = (a^sp'*',0,ksr) , 

p' = (0,p'“,0) = (O,pp",0) , 

K = (0: x[p'~, k'lr) , ^2 = (0, X 2 p'~, k'2T) , fcg = (0, , k'gr) , 

q = {q^,q~,o) = (p+,r^p",o), 

qi = {yq^,yq~, 92 = ((i-p)g^,(i-?/)g“,-qT), ( 5 ) 


with the constants, 


p = I - , 


Ma, ■ 


( 6 ) 


ki {k[) is the b (s) quark momentum, Xi (x'P) are the momentum fractions associated with the Ab (A) baryon, and 
k^^ the corresponding transverse momenta, satisfying k|^ =0. y is the momentum fraction and qr the transverse 
momenta carried by the charm quark in the J/ip meson. 

The structure of the A;, baryon wave function Ta;, is simplified under the assumptions that the spin and orbital 
degrees of freedom of the light quark system are decoupled, and that the Ab baryon is in the ground state (s-wave). 
The wave function is given, in the transverse momentum space, by 




2V2N^ 

fA, 

8V2Nc 


n 

1^2 




(27r)3 


e*"‘-“"e“'’"(0|r[6“(0)4(u;2)d;(u;3)]|Ab(p)) , 


[{^ + MAth5C]f3j[Ab{p)]a'i'{ki, p) , 


(7) 
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where b, it, and d are the quark fields, a, b, and c the color indices, a, (3, and 7 the spinor indices, C the charge 
conjugation matrix, Ab{p) the A;, baryon spinor, and /a^ the normalization constant. 

The A baryon wave functions are defined, in the transverse momentum space, via 


{YAU^{k',p) = I , 

-^■^^^(o'M>^p'‘'C')/37[7''75A(p')]a$'^(fc-,Ai) , (8) 

with the normalization constants /a and /J, the A baryon spinor A(p'), and the definition = [ 7 ^, 7 j/]/ 2 . The 
wave functions and are associated with the different spin structures of the three valence quarks in the A 

baryon. 

The J/ip meson wave function is expressed as 

1 f d^in 

{YJ/^U{q^,X,p) = ^e*^--(0|T[cMO)c„H]|J/V^(g,A)) 

= +Mj/^)]a/3U{qi,fJ,) , (9) 

where the J/4> meson decay constant /j/^ has been absorbed into the wave function 11 , and £^( 9 , A) is the polarization 
vector of the J/'ip meson with the helicity A. It has been assumed that the J/ij} meson wave functions associated with 
longitudinal and transverse polarizations possess the same form in the above expression. To ensure that the valence 
charm quarks are close to the mass shell, 11 should have a sharp peak at ?/ ^ 1/2. The J/i/ meson wave function 
has been determined from the data of the B decays 0. For the factorizable amplitudes, the momentum 

fraction y is integrated out, and only the decay constant /j/^ is relevant. 

When the transverse degrees of freedom of partons are taken into account, the factorization of a QCD process 
should be constructed in the b space with b being the variable conjugate to kr- This is the reason 1/b serves as 
the factorization scale stated in Sec. II. Sudakov resummation of the double logarithms and RG summations of the 
single logarithms contained in the above hadron wave functions lead to 


= exp 

,bi,y) = exp 
^{qt^bg,y) = exp 


J2s{w,k^)-3 f ^j{as{y)) 
1=2 ^ 

2 

- 3 / ^'y{as{fi)) 
J w' M 


Z =1 

2 


'^s{wq,q+)-2 f ^7(a^(/r)) 
1=1 “'“<1 ^ 




T^{y), 


( 10 ) 


with the superscript j = V, A and T, and the quark anomalous dimension 7 = —aaj'K. The Sudakov exponent s for 
the J/ijj meson depends on the dominant component g/". The factorization scales w, w' and Wq are chosen as 


with the variables, 


w = min 




w = mm 


1 1 1 



bi = |b2 - bal , 


6 ; = |b' 2 -b' 3 |. 


( 11 ) 


( 12 ) 


The initial conditions ip, (j)^ and tt of the Sudakov evolution absorb nonperturbative dynamics below the factorization 
scales w, w' and Wq, respectively. Note that the Sudakov effect from the J/ip meson is weak, since the distribution 
amplitude TT{y, A) vanishes rapidly as ?/ —> 0,1. 

The explicit expression of the Sudakov exponent s is given by 


s{w,Q) 


— ln—A{asip)) + B{as{p)) 
PIP 


(13) 
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where the anomalous dimensions A to two loops and B to one loop are 


A = Cf—- 

TT 

2 ttg 

B = -^ In 

3 TT 


[67 

TT^ 

10 

8.3 n eA ®1 


. 9 

3 

27 " ' ' 

—Po In- 

3 2 

(v) 


g27E-l 


(14) 


Cf = I being a color factor, n/ = 4 the flavor number, and 7 f the Euler constant. The one-loop running coupling 
constant. 


aAf^) _ I _ 

TT ^oHm^/Aqcd) ’ 


with the coefficient /So = (33 — 2n/)/12, will be substituted into Eq. (13). 


(15) 


IV. FACTORIZATION FORMULAS 


The Af, ^ J/ipA decay amplitude A4 is expressed as 

M=i^M{\) , (16) 

A 

where the amplitude A4(A) can be decomposed into different structures with the corresponding coefficients Ai, A 2 , 
Bi, and B 2 '. 


Ai ^75 + ^ 2-^—75 -I- Bi^ + 


A&(p) . 


M{X)=A{p') 

The general factorization formula for A4(A) is written as 

MiX) = J[Bx] J 

X (^J/i/-)pp' («> 2/> bg, A, fl) Ma,,p) iYA,)a0jip, X, 5, n) , 

with the measure of the momentum fractions, 

[Bx] = [dx][dx']dy , [dx] = dxidx2dx^5 [ 1 — ] , [dx'] = dx'idx^dx'^S j 1 — j . 


(17) 


(18) 


(19) 


Z=1 


Z=1 


The measure of the transverse extents \ph] will be defined in Appendix A. 

The RG analysis of H leads to 

H{x,x',y,b,b',bg,MA^,p.) = exp 

xn{x,x',y,b,b',bq,MAt,t) , (20) 

with the integer n = 6 for the factorizable diagrams and n = 8 for the nonfactorizable diagrams. The superscripts a', 
P', • • •, have been suppressed. The argument t in the initial condition of H implies that the Wilson coefficient and the 
running coupling constants as for the hard gluons are evaluated at t. This initial condition will be calculated based 
on the lowest-order diagrams in Fig. 1. For the purpose of presentation, we shall rewrite the initial hard amplitude as 
Ti = aHpQ., where a is the Wilson coefficient, Hp the numerator of Ti. depending on the spin structure of the three 
valence quarks in the A baryon, and the Fourier transformation of the denominator of Ti. from the kp space to the 
b space. 

Substituting Eqs. and (H^) into Eq. (|l^, we derive the factorization formulas, 

= 2GpV,,V:, Y. / \DbX\^s{f)fa\f) 

X'ij:{x)(IP {x')'K{y) exp[—SXF[pXV‘ , ( 21 ) 


—n 


dp 

a M 
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where F represents the coefficients A 2 , Bi, and B 2 , the superscript i labels the diagrams in Fig. 1, and the 
superscript j labels V, and T associated with the spin structures of the valence quarks in the A baryon. The factor 
2 comes from the symmetry under the u-d exchange. The Wilson coefficients are a* = Ci/Nc for the nonfactorizable 
diagrams with only one hard gluon attaching the c quark or the c quark, and a* = 02 = C 2 + Ci/iVc for the factorizable 
diagrams and for the nonfactorizable diagrams with two hard gluon attaching the two quarks of the J/tp meson. At 


the characteristic scale of order Ci/W is about few times larger than a 2 - The Fourier 

integrals employed in this work, and the resultant measures \Db\' of the transverse extents for each diagram i are 
given in Appendix A. The explicit expressions of are presented in Appendix B. The functions F[^p are given in 
Appendix C. To derive the expressions of i?];?, we have chosen xi = 1 and X 2 = X 3 = 0, since X 2 and X 3 , being of 
0 {A/mb), are negligible in the current leading-power analysis. 

The main theoretical uncertainty of our predictions come from higher-order corrections to the hard amplitudes, 
which can be reflected by the choice of the hard scale f We shall consider the following two choices: 


f = max(t 5 ^, ^ 2 , w, w', Wq) , 


f = max 


-,W,W ,Wq 


( 22 ) 


where the hard scales t\ and ^2 associated with the two hard gluons in each diagram i are listed in Table I. The 
max in the above expressions simply means that the hard scales should be larger than the factorization scales. It is 
expected that the first choice will lead to values lower than those from the second one. The theoretical uncertainty 
can be reduced after next-to-leading-order corrections are included. 

The exponents S'* are given by 


S* = ^s(w,x/p+)-f 3 / ^7(as(/r)) 


1^2 


+ ^s{w\x[p' )-f3/ ^-f{as{p)) 


1^1 


S* = ^ s{w, xip^) + 3 I ^-/(asin)) 

1=2 ^ 

+ '^s{w',x'ip'~) + 3 f ^7(q!s(/I)) 
J w' M 


1^1 

2 


+ ^s{wq,yiq+) + 2 I ^j{asip-)) , 
1=1 '^'^1 ^ 


(23) 


(24) 


for the factorizable and nonfactorizable diagrams, respectively. The Sudakov exponential associated the J/ip meson 
distribution amplitude appear only in the nonfactorizable diagrams. The evolution factors described by 7 correspond 
to g{t,b) introduced in Sec. II. 


V. NUMERICAL RESULTS 


In this section we evaluate the factorization formulas in Eq. 
amplitude tp, we adopt the model proposed in p^], 


(Eh) numerically. For the A^ baryon distribution 


p){xi,X2,X3) = NxiX2X3exp 


_ rnf 

2f3"^xi 2/3^x2 2 / 32 x 3 


(25) 


where the shape parameter /3 = 1.0 GeV, and the mass of fight degrees of freedom in the baryon, mi = 0.3 GeV, 
have been determined in Q . The normalization 


J [fia:]^/i(a:i,a;2,a;3) = 1 , 


(26) 
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leads to the constant N = 6.67 x 10^^. The constant /a^ = 2.71 x 10“^ GeV^ and the above Af, baryon distribution 
amplitude tp have been chosen to satisfy the experimental upper bound of the A;, —> AJP branching ratio and the 
heavy quark symmetry |^. 

The A baryon distribution amplitudes have been derived using QCD sum rules |ll| , |l2|. The asymmetric distribution 
in the momentum fractions of the three quarks implies SU{3) symmetry breaking |12(| . In this work we adopt the 
models proposed in 0 , 

(j)^{xi,X2,X3) = 42(j)asixi,X2,X3) [ 0 . 111 (a ;3 - xl) + 0.093(a;2 - xa)] , 

(j)^{xi,X2,X3) = -42(/)as(a;i,a;2,a;3) [0.093(a;2 + xl) + 0.376x1 
—0.194x2X3 — 0.207xi(x2 + X 3 )] , 

(/>^(xi,X 2 ,X 3 ) = 42^as(xi,X2,X3) [3.6 (x2 - X 3 ) + 0.32(x2 - X 3 )] , (27) 

with the asymptotic distribution amplitude, 

()ias(Xi,X 2 ,X 3 ) = 120 X 1 X 2 X 3 . (28) 

The normalization constants /a and are chosen as 

/a = 0.45 X 10“^ GeV^ , /J = 0.2 x 10"^ GeV^ . (29) 

It is easy to observe that the above A baryon distribution amplitudes satisfy the relations, 

()i^(xi,X2,X3) = -0'^(Xi,X3,X2) , 

()|^(X1,X2,X3) = ()>^(X1,X3,X2) , 

()|'^(X1,X2,X3) = -^'^(X1,X3,X2) , (30) 

and the normalization, 

J[dx](l)^’^’'^{xi,X2,X3) = 1 . (31) 

It is most likely that the two charm quarks in the J/ip meson carry the equal fractional momenta, such that they are 
close to the mass shell. Hence, the following J/ip meson distribution amplitudes with sharp peaks at the momentum 
fraction y = \ 12 have been proposed |^], 


= — fj/^ 2/^(1 - yf , (32) 

with the J/pj meson decay constant fj/^ = 390 MeV [^. It has been found that Eq. ( |^ ) gives the branching ratios 
H(i? —> J/ipK^*'^) in agreement with experimental data 1^. We have confirmed that the PQCD predictions are not 
sensitive to the functional form of tt, as long as it has a sharp peak at y ^ 1/2. 

Employing the parameters Gp = 1-16 x 10“® GeV“^, Vcb = 0.04 and Vcs = 0.975, Aqqd = 0.2 GeV, and the masses 
= 5.624 GeV and Mj/^ = 3.097 GeV, we obtain the coefficients Ai, A 2 , Bi, and B 2 . The contributions from the 
transverse component of the A baryon distribution amplitudes are listed in Table II. These coefficients, dominated 
by the nonfactorizable amplitudes, are mainly imaginary. The Af, ^ AJ/ip decay rate is written as 


r = 


1 Pc 
SttMI 


\M\\ 


(33) 


where Pc is the magnitude of either of the final-state particle momentum in the center-of-mass frame. Following the 
choices of the hard scale t in Eq. (^^, we derive the branching ratio, 

B{Ab AJ/Pj) = (1.7 - 5.3) X 10"^ , (34) 


for the Af, baryon lifetime r = (1.24 ± 0.08) x 10 s. Here we do not consider the minor uncertainty from r. The 
value in Eq. (|^ ) is consistent with the experimental data [p7| , 

B{Ab AJ/iP) = (4.7 ± 2.8) x 10"^ . (35) 
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We also calculate the asymmetry parameter a associated with the anisotropic angular distribution of the A baryons 
emitted in polarized Ab baryon decays: 


dr oc (1 + ap' ■ V) , 

V being the A;, baryon polarization. The explicit expression of a is given by |pl 

+ 2E]i^Re[{S + DyP,] 


with the factors, 


2(|^P + \P2?)M^j,^ + (|5 + D\^ + 


(36) 


(37) 


S = -Ai, 
D = - 


Ej/ip{E]\ + Mi\) 


(^1 — A2 ) , 


P = - 


Pc 


E 


P 2 = 


j/ip 

Pc 


Ea + Ma 


AfAfc + ATa 
Ea + Ma 

Bi , 


Bi + B 2 


(38) 


where Ej/^ (Ba) is the energy of the J/y meson (A baryon), and the A baryon mass Ma has been set to zero for 
consistency. The result is, following Eq. (22), 


a =-0.17 ~-0.14 , (39) 

consistent with those derived from the quark model based on FA shown in Table III. The predictions for a are 

very stable with respect to the variation of the hard scale t and of the hadron distribution amplitudes. Therefore, it 
serves as an ideal quantity to test the PQCD approach. 


VI. CONCLUSION 

In this paper we have analyzed the nonleptonic heavy baryon decay A{, —> AJ/ip using the PQCD formalism. It 
is a special heavy-to-heavy mode, to which PQCD is applicable, because of the soft cancellation between a pair of 
nonfactorizable diagrams. We have shown that nonfactorizable contributions dominate in this mode, similar to the 
heavy meson cases B J/yK^*\ since factorizable contributions are suppressed by the small Wilson coefficient 02 . 
After considering the theoretical uncertainty arising from higher-order corrections, we have derived the branching 
ratio B{Ai, —> AJ/ip) and the asymmetry parameter a associated with the anisotropic angular distribution of the 
A baryons produced in the polarized A^ baryon decays. The latter quantity is stable with respect to the variation 
of the hadron distribution amplitudes and to higher-order corrections. The comparison of our predictions with the 
experimental data will provide a test of the PQCD formalism. 

The PQCD results are consistent with those derived from other approaches based on FA as indicated in Table 
III. However, we emphasize that the theoretical bases between PQCD and FA are very different. In PQCD, 02 is 
a Wilson coefficient, which is small in a wide range of the energy scale. To explain the experimental data, large 
nonfactorizable contributions are necessary. In FA nonfactorizable contributions are neglected. To account for the 
data, 02 must be treated as a free parameter, and the value 02 ~ 0.23 0 has been adopted. Also, nonfactorizable 
contributions are imaginary and their strong phases can be evaluated in PQCD, while the parameter 02 is real (or 
with an arbitrary phase) in FA. As stressed in the imaginary nonfactorizable amplitudes determine the relative 
phases of the various B Dn decay modes, which are essential for explaining the recent Bd —> data pl|,p^. 


This work was supported by the National Science Council of R.O.C. under Grant Nos. NSC-90-2112-M-001-077, 
NSC-90-281 l-M-001-042 and NSC-90-2112-M-001-038. 










APPENDIX A: FOURIER INTEGRATIONS AND B MEASURES 


We list below the Fourier integration formulas that have been employed in the derivation of the hard amplitudes 
in the b space, where Ji, IVi, Kq and Ki are the Bessel functions, and Zi the Feynman parameters: 


—7 = 2'KKo{\/Ab) , A > 0 , 


P + A 


(Pk 


„zkb 


Jo 


(fc 2 + A)(F + i3) " Jq 

gi(ki-bi+k 2 'b 2 ) 


5 ■*•5 


A, B > 0 


d^kid^k2 


{kl + A){kl + B)[{ki + k^f + C] 


dzidz2 


Mix, 


(vW^)d(Z 2 ) 


TT 

+ 2 


0 zi{l-zi) ^/\^\[ 

Nl{\/X2\Z2\) — iJl{\/X2\Z2\) d( —2'2)|' , 

A > 0 , and B, C arbitrary , 

i(kibi+k2b2+k3b3) 


nl 


(fc2 + A)(fc| + B){kl + C)[{kx +k2 + kzY + D] 
Jo zi(l -Zl)z2(l - 22 ) [ 


TT 

2 I 


Nii^/XolZol) - ei-Zo) 

A, B > 0 , and C, D arbitrary , 


with the variables, 


Zi = A z -\- B {1 — z) , 
Z 2 = A(1 — 22 ) + 


2:2 


zi(l - Zi) 


[B (1 — zi) + C Zi] 


V fu k ^2 , ^ 1(1 “ ^l)k2 

X2 — (oi — Z1O2) H--O2 , 


2:2 


■^3 — A (1 — Z 3 ) + 


23 


22(1 - Z 2 ) 


B {1- Z 2 ) + 
Z 2 O - Z 2 ) , 


2^2 


Zi(l - Zi) 


[C (1 — zi) + D zi] 


X 3 = [bi - b2Z2 - &3-Z2(l - Zi)f H- -{b2 - bozif 


Z 3 


Zi(l - Zi)z2(l - Z2) 


bi. 


Z2Z3 


After performing the above integrations, we combine the 2 : measures with the ordinary b measures, 

dHi 


[dbi] = 


{2ny 


to form the special b measures [Db] appearing in Eq. (©• [Bby for each diagram i in Fig. 1 is given by 

= [Vb]^d) = {2TTf[db2][db'^][db'o]dzidz2 , 


(Al) 

(A2) 


(A3) 


(A4) 

(AS) 

(A6) 


(A7) 


(AS) 


(A9) 


[Vbp = [Vbp = = [Vb]^'^ = = [Vb]^^'^ = {2nf[db2][db'y^[db'o][dbq] , (AlO) 
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= [2'K)^[dh2][dh'^[dh'^dzi , 


(All) 


= {2T:f[dh’2\[d%][dhq]dzidz2 , (A12) 

= {2TTY[db'2\[db''^[dh^dzidz2dz^ , (A13) 

p6]('=) = = [l?6](’-) = {2TTf[db2][db3][db'^][db'^] . (A14) 


APPENDIX B: EXPRESSIONS OF 0.^ 

Employing the integration formulas in Appendix A, we derive the hard amplitudes in the b space. The hard function 
12* for each diagram i in Fig. 1, which arises from Fourier transformation of the denominators of the internal particle 
propagators, is expressed as 

(y/^) , (Bl) 


with 

Ba = - 6')2 + [(1 _ ^ 

2^2 

Za = X2X2pMl^{l - Z2) + ^^^^^^^^^ [(1 - a;i)2:i + 2:3(1 - X2)(l - Zi)]pMl^ , 

(B2) 


with 




Ko [yX2x'2pMK^{b2 + bq)^ Ko (\/a;3X3pMA^&3^ 

xKq ^^Jx3{l - X2)pMA^\b2 - 62 !^ [y^bMf^^b^ 6{Zi,) 



r‘2 

Zb = — -{ 1 -X 2 - y)[xi 2 P+{l - y)r‘^] , 


(B3) 


(B4) 


with 


12W = Ko (^\/X2X2pMAt\b2 -bg\J Kq (yxox'^pMAi,b'o 

xKo (^^Jx3{l-x'2)pMA,\b2 - b' 2 \^ j^o (v^Ma,6,) e{Z,) 
-Nq (^^/\Zc\MAhbq^ + iJo (^VJz^MA^bq^ 0(—Zc)! , 


TT 

'2 L 


Zc = — + ix2 - y)[x2P + (1 - y)r‘^] , 


= 2 ^^ ^° {yX 2 x' 2 pMA^b 2 ^ Ko [yXox'opMA^b'^ Ki (^a/^^aJ &2 - ^ 2 !) > 


with 


(B5) 


(B6) 

(B7) 
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with 


with 


with 


Zd = [a ;3 + 2:2(1 - X 3 )p\MI^Zi + 0:3(1 - x' 2 )pM‘l^{l - zi) 


+ K)) \k, ( x /^) 9{Z,) 
+ 1 [iVi (^V\Ze\Be) - iJl (Vl^el^e)] ^(-^e)| , 


Zi(l - Zi) 2 


+ (^2 + ? 


Ze = X 2 X 2 pMj^^{l - Z2) -\- ^— - {xi - y){l - x[p - yr )J Zi 

+ ^ - X3 - y){x3P + (1 - y^) ( 1 -zi)|m^^, 

^Ki (Jz7b 7) e{Zf) 
8 zi(i-zi)z 2 (i-^ 2 ) f f) ^ 

+ 1 k U\Zf\B^-iJ, 0(-Z/)| , 


^ Z3(1-Z3)Z2(1-Z2) ^2 ^ + (1 _ 


+ ^'2 - bqZl - [b's + {I - Zl)hq]z-^ , 


Z/ = X 2 X 2 pMl^{l - Z 3 ) + ^^(l_^^) k3^3P^A,(l - Z 2 ) 


Z2 \ r 


Zi(l - Zi) 


- (1 - 2:3 - yi)ixsP + (1 - yi)r^)^ (1 - Zi) 


+ i — + {x2-yi){x2P + yir'^)jzi I , 


[\/X2x'2pMp,^b‘^ Kq k^3a:3pMA^(&3 + 5,)^ 

X-ftio ^^3 + ^ 2(1 - X3)pMA^\b2 - 62 !^ ^(^9 


+ - -NoU\Zq\MA,bq\+tJoU\Zq\MA,bq\ 0(-Zg)L 


^9 = ^ - (1 - 23 - J/i)[a:3P + (1 - yi)r‘^] 


(a/^^MaJ63 - 6 , 1 ) (a/^) 0(^^) 

+1 k (v/|^,,lk -zJi (vW^)] e{-Zh) \, 


(B8) 

(B9) 

(BIO) 

(Bll) 

(B12) 

(B13) 

(B14) 

(B15) 
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with 


with 


with 


Bu = ^ ^ 


22 
„' „ A/f 2 


Zh = X2X'2PMI^{1 - Z 2 ) + 


Z2 


Zi{l - Zi) 


+ (1 - XI - yi)(l - x[p - (1 - 2 /i)r 2 ) 


21 


-^ + {x3 - yi){x3P + yir'^) 


{l-z,)\Ml , 


= Ko [\JX2x'2pMi^^b'^ Ko (\/x3X3pMaJ&3 -bq\^ 

xKo (^^X 3 +x^( 1 -X 3 )pMaJ &2 - & 2 I) 1^0 (\/^Ma, 6 ,) 
-TVo +zJo 0(-^*)| , 


7 T 

2 L 


= ^ + (a;3 - y)(4P + yr"^) > 


= cnrriT) (^/ 


= ^i^^-^(62 - 6^)^ + [6^ - 21(62 - 5')]2 , 


ZjBj 1 , 


22 

7- = 

' 2i(l-Zi) 


|(1 - x'i)(l - Zl)p + [X3 + X2(l - X3 )p]2i|m^^ 


+X 3 X 3 PM? (1 - Z 2 ) 


fl(^) = Ko (1 - xi)(l - x[)pMjy^h'^ Ko (\/x3X3pMaJ62 - 63 ]) 

xKo (^\J (1 - xi)x 3 pMaJ 62 - 62 - 63 + 63 !^ ifo (\/(1 - a^i)p^Aj 63 - 63 I) , 

= Ko (1 - xi)(l - x()pMaJ62 - 63 + 63!^ Ko [yX 3 x'opM^^\b '2 - 63I) 

Ko (^\Jx 3 {l - X3)pMaJ62 - 62 - 63 + 63!^ Ko - Xi)pMA^\b 3 - 63I) , 

= i^o (^\/ (1 - a;i)(l - x[)pMA^{b2 + bg)^ Ko (\/x^^MaJ62 - 63I) 
xKo (^yi - xi)x'opMA,\b2 - 6'2l) 1^0 (\/^Ma,6,) e{Z^) 

-Ko (^\/\Zm\MAi,bq^ -\-iJo (^\/\Zm\MAi,bq^ d{ — Zjn) 


TT 

'2 L 


(B16) 


(B17) 

(B18) 

(B19) 


(B20) 


(B21) 


(B22) 


(B23) 


with 
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(B24) 


with 


with 


with 


j,2 

Zm = ^ + {yi - a:i)[l - j/ir^ - x'^p] , 

= Kq (^sJ (1 - a;i)(l - x[)pMA^{b2 + Kq (\/^3^MaJ62 - &3I) 
xKq ^ya;3(l - x[)pMAt\b2 - 62 !^ |^o 6 »(Z„) 

+ 1 [-A'o +*^o (v 1 ^Ma, 6,)] 0(-Z„)| , 

^2 

^« = ^ + (j/i - a;i)[l - yir^ - a;'ip] , 

= Ko (1 “ ^i)(l “ Xi)pMAt\b2 - 5,1^ Ko (\/x^^MaJ62 - 63I) 
xKo (^^(1-xi)x^pMaJ62 - 62I) 1^0 (v^Ma,6,) 

+ '^ “-^0 ^\/|-^o|-/WA6&g^ + * "/o ^\/|-^o|AfAjfog^ 0( —Zo)|’ , 

Zo = ^ + {I - Xi - yi)[l - (1 - yi)r'^ - x[p\ , 

n^p') = Ko (1 - Xi)(l - x[)pMA^\b2 - bqlj Ko J^2 “ ^sl) 

Ko (^^Jxo(l-x'^)pM^,\b2 - b'^\^ |xo (v^Ma,6,) e{Zp) 

^0 ^\J I ■^P I -^Ab + * •A) ^ I •^PI -^Ai, bq^ 9[ -^p)^ ) 

p2 

= ^ + (1 - a:i - 2/i)[l - (1 - 2/i)r^ - x[p] , 

17(9) = JsTo (1 - xi)(l - x[)pMA^\b2 + 63 - Ko (\/^E3^^aJ& 2 - feal) 

XitTo (1 - Xi)x3PMaJ62 - &2 - ^3 + ^3!^ ^0 (a/(1 - Xi)pMA^\b3 - 63!^ , 



(B25) 

(B26) 

(B27) 

(B28) 

(B29) 

(B30) 

(B31) 

(B32) 
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APPENDIX C: EXPRESSIONS OF Hy 


The part H^p of the hard amplitudes, that varies among the spin structures of the three valence quarks in the A 
baryon and among the form factors Ai 2 and Bi 2 , are gathered in the following tables. Their derivation is similar to 
that presented in for the proton form factor. The contributions to from Figs. 1(a), l(k), 1(1), l(n), l(p), and 
l(r) vanish. 


F 

3 

hWj 


V 

-8r(l-r^)^MX,(l-4) 

A 

8 r(l - rrMyll - 4) 


T 

24r(l-r^)^MX,(l-a:^) 

A 2 

V 

16r(l - r^)My{l - - 2y) 

A 

-16r(l - - x^il - 2y) 


T 

-48r{l-ryMll{l-xy{l-2y) 


V 

8r(l-r^)-MX,(l-4) 

Bi 

A 

-8r(l-ry^Ml{l-xy 


T 

—24r(l — FyMy{l — X 2 ) 


V 

16r(l - r^Myil - x^il - 2y) 

B 2 

A 

-lQT{l-ryMiyi-x'2){l-2y) 


T 

-48r(l-r")MX;(l-x')(l-2y) 


F 

3 


Ai 

V 

16r(l - r^)^MXJl - 4 ) 2 / 

A 

-16r(l - r^y^Miyi - x^y 


T 

-48r(l - - x^y 

A 2 

V 

0 

A 

0 


T 

0 


V 

-16r(l-r^)^M4(l-4)2/ 

Bi 

A 

16r(l-r^)^M4(l-4)2/ 


T 

48r(l - r^)^M4(l - x^y 

B 2 

V 

0 

A 

0 


T 

0 


F 

3 

jjWj 

Ai 

V 

16r(l — — X 2 ) 

A 

— 16r(l — r^)^M4(l — X 2 ) 


T 

-48r(l-r^)^M4(l-4) 


V 

0 

A 2 

A 

0 


T 

0 


V 

— 16r(l — r^y^Myy — X 2 ) 

Bi 

A 

16r(l — r^)^M^^(l — a;^) 


T 

48r(l-H)^M4(l-4) 

B 2 

V 

0 

A 

0 


T 

0 
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4r(l - (6 + - 2x'^ + - 2^^ 

— 12y — 4r^y + ‘ix'^y — Ar'^x'iy + 4y^ + Ar'^y'^) 


—8r(l — r (5 + + X3 — r^a;^ — lOy 

— 2r^y + 2x[y — 2r‘^x[y + 4y^ + 2r^y^) 



—4r(l — r"^)M'^^{—2 + — 2xi + 2r"^Xi + 2x^ — 2r^x^ 

— Ar'^y + 4y + 4:x[y — Ar'^x'^y — 4y^ + 4r^y^) 


—8r(l — r (3 — r + Xi — r'^x^ — x^+ r x'^ — 6y 

+ 2r^y — 2x[y + 2r‘^x[y + 4y^ — 2r^y^) 


16rM'^^{l — 3xi + 6r^x\ + Xg — r^Xg — 2r^y + 4xgy 
— 4r^Xgy — 2 y^ + 2 r^y^) 


— — 3xg + 3r^Xg + Xg — r Xg — 2r^y + 4xgy 

— 4r^Xgy — 2y^ + 2r^y^) 


16rAf^^ (1 + r + 3xg — 3r Xg — Xg + r^Xg — 2r^y — 4y 
— 4xgy + 4r^Xgy + 4y^) 



8r(l — [3 — 5y + 2r y + 2x^y — 2r^x^y + 2y^ — 2r‘‘y 


- 32 r{l-r^)Ml{l-yr 


1 - 2 / 


16r(l-r^)M? (3-2y)(l-y) 


-8r(l — r^)M^ (—3 + by + 2r‘^y + 2x^y — 2r x^y — 2y^ — 2r^y^) 


32rMl{l - y)(l -y- r^y) 


32rMX, (1 - y)(-l + y + r^y. 


l&rMi (—3 + 5y + 5r y + 4x3y — Ar x^y — 2y — 6r y ) 


-8r(l - r^)Ml (2 + X2 - r^x^ - 2y - 2x2y + ^r^x'^y) 


8r(l — r )M? (4 + 2r^ + X2 — r X2 + 2x3 — 2r Xg — 4y — 2r y — 2x2y + 2r^X2y) 


_ 16r(l(-2+ 2 y + x^)_ 


16r(l-r^)A4X (4-x^-4y) 


-8r(l - r^)Ml(2 - X2 + r^x'^ - 2y + 2x2y - 2r^X2y) 


8r(l - r"^)MX, (2 - X2 + r"^X2 - 2y + 2x2y - ^r^x’^y) 


-8r(l — r )M? (—4 + 2r + X2 — r x'^ + 2x^ — 2r‘^x^ + 4y — 2r^y — 2x2y + 2r"^X2y 


16rMX (2 — 3x2 + 3x^X2 — 2y — 2r y + 4x2y — 4r^X2y) 


— 16rMX (2 — 3x2 + 3r X2 — 2y — 2r y + Ax^y — 4r^X2y) 


16rMX (—4 + 2r^ + 3x2 “ 3x^X2 + 4x3 — 4x Xg + 8y — 4x y — 4x2y + 4x^X2y) 
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Bill 


Ai 




— 8 r(l — r (3 + 3r^ — 4 x 3 + 4 r^X 3 — 2y — Ar^y) 

A 2 

H 

16r(l — r^)M^^(l — 2 x 2 y) 

S 

-16r(l - r^)Ml,(l - 2 a:^y) 


H 

16r(l - r^)M£ (-3 + 2 y + 2x^y) 

Bi 

H 

8 r(l — r"^)M^^(l — + 2 x 3 — 2 r^X 3 + 2r^y) 

El 

— 8 r(l — r^)M^^(l — + 2 x 3 — 2 r^X 3 + 2 r^y) 


m 

— 8 r(l — r‘^)M\^ (3 — 3r^ + 4^3 — 4 r^X 3 — 2y + ir'^y) 

B 2 

m 

-16r(l-r^)MX,(l-2a;^y) 

m 

m 

16r(l — r^)M^^(l — 2 x 2 y) 


16r(l-r^)M^^{5-2y-2x'2y} 
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Table I: The hard scales t\ and for each diagram i in Fig. 1. 


i t\ _ 

(a) max( y/x2x'^pMK^, ) 

(b) max( y /x2x'^p MA„ ) 

(c) max( ^X2x'^pMa„ ) 

(d) max( ^X2x'2pMa,, ) 

(e) max( ^/x^j^Ma,, ) 

(f) max( y/x2x'2pMA,, ) 

(g) max( ^X2x'2pMa„ ) 

(h) max( y/x2x'2pMA,, ) 

(i) max( y/x2x'2pMA,, ) 

(j) max( y^X2x'2pMA,, ) 

(k) max( a /(1 - a;i)(l - 

( l ) max( a /(1 - a;i)(l - x[)pMa^, 

(m) max( a /(1 - a:i)(l - x[)pMAt, 

(n) max( a /(1 - a;i)(l - x'J/jMa,, 

(o) max( a /(1 - a;i)(l - x'J/oMaj, 

(p) max( a /(1 - a;i)(l - x[)pMa^, 

(q) max( a /(1 - a:i)(l - x'JpMa,, 

(r) max( \/(l - a;i)(l - x[)pMa,, 


_^2_^_ 

max( y/xsx'^pMA,, ) 

max( i/xaXgpMAt, ) 

max( y^xsx'^pMA,, ) 

max( i/xaXgpMAt, ) 

max( y ^xsx'^p MA,, ) 

max( ) 

max( a/xs^pMa,, 
max( yJx^x'^pMAt,, 1^/^^ ) 
max( y/xsx'^pMAt, ) 

max( y^xsx'spMA^, ) 

) max( a /xs^p Ma,, ) 

|b,-bUb^| ) ^ /x3x'3p MA„ ) 

) max( ^xsx'sp Ma,, ) 

) max( y ^xsx'sp MA,, ^ ) 

) max( ^xsx'sp Ma,, je^ ) 

) max( ^x3x'^p Ma,, ^ ) 

| b '- 4+63 1 ) Vxsx'sP Ma,, je^ ) 

]4t) V^sx'spMa,, i^T ) 
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Table II: The coefficients Ai, A 2 , Bi, and B 2 from the transverse component (jF of the A baryon distribution 
amplitudes for the choices of the hard scales, 

(a) t* = max (t^, w, w', Wq), (b) = max , w, w', . 



(a) 

(b) 

A{ 

-7.17 X 10-1^-8.69 X lO-^i 

-1.18 X 10-y- 1.54 X 10-«i 

^2 

-5.91 X - 1.62 X 

5.97 X 10-1U-2.83X 10-«i 

Bi' 

6.54 X 10-1^ - 7.53 x IQ-^i 

1.14 X 10-y - 1.27 X 10-«i 

Bi' 

2.63 X 10-11 + 1.21 X 10-«i 

6.36 X 10-1^ + 1.83 X 10-«i 


Table III: The Af, - 
data are also shown. 


AJ/il) branching ratios and the asymmetry parameters a derived in different approaches. The 



Ref. PI] 

Ref. 1^ 

Ref. |2^ 

This work 

Experimental data 

B(Ab ^ AJ/^) 

2.1 X iir^ 

6.04 X K)-4 

2.49 X 

(1.65 ~ 5.27)x 10-^ 

(4.7 ±2.8) X 10-^ 

a 

-0.11 

-0.18 

-0.208 

-0.17- -0.14 

- 
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FIG. 1 

































